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PREFACE 


A previous report. Reference 1. prepared under this contract, described 
a computer program, designated RETSCP. for the analysis of rocket 
engine thrust chambers with cyclic plasticity. That report included 
a complete description of the program logic . together with a users 
manual. It is the purpose of this report to illustrate the detailed 
application of the RETSCP program. Volume I describes the analysis 
of a regeneratively cooled OFHC copper combustion chamber. 

Operating conditions are such that plastic strains dominate and 
fatigue life is in the low cycle regime. Volume II describes the analysis 
of an attitude control thrust chamber. Operating conditions for that 
engine are such tha^ the material behavior is elastic . or on the 
threshold of the elasto-plastic regime; thus . the fatigue life is in 


the high cycle regime 



SUMMARY 


A finite element stress analysis was performed for the film cooled 
throat section of an attitude control thruster. The analysis employed 
the RETSCP finite element computer program (NASA CR-134640). The 
analysis included thermal and pressure loads, and the effects of 
temperature dependent material properties , to determine the strain 
range corresponding to the thruster operating cycle . The configuration 
and operating co.idltions considered, correspond to a flightweight 
integrated thruster assembly which was thrust pulse tested at the 
NASA Lewis Research Center. The computed strain range was used 
in conjunction with Haynes 188 Universal Slopes minimum life data 
to predict throat section fatigue life . The computed number of cycles 
to failure was greater than the number of pulses to which the thruster 
was experimentally subjected without failure. 
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INTRODUCTION 

As part of the NASA goal to use all cryogenic fueled engines for the re- 
usable Space Shuttle « the necessary technology for reliable attitude 
control thrusters must be developed. These thrusters are unique In 
that the life requirements are quite demanding. These 1500 pound 
(6672 N) thrust engines « using low temperature oxygen and low temperature 
hydrogen gases, must be capable of firing (pulsing) a minimum of 500 
times per mission for 100 missions. As part of this technology effort 
It becomes Imperative that an analytical design procedure be formulated 
to Insure that the fatigue life goal Is met under thermal and pressure 
cycling . 

Combustion chambers can be designed to operate at suitable temperature 
levels by employing film cooling. The thermal analysis of cooled 
combustion chambers follows established methodology. Empirical data 
Is available for estimating convection heat transfer. The temperature 
distribution within the structure Is then computed using one of many 
available finite difference computer codes for transient or steady state 
heat conduction In two or three dimensions . 





The state of stress In film cooled rocket chambers varies In three 
dimensions . The situation Is further complicated by the dependence of 
material properties on temperature. It Is the range of strain, through 
which a given chamber Is cycled, that determines fatigue life. 
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In order to gain experimental data for such attitude control thruster 
configurations, accelerated testing Is attractive. In this manner, 
fatigue testing Is carried out at conditions which are more severe 
than the normal operating conditions. The resulting stress field may 
be elastic or may be elasto-plastlc . 

In this report, a method Is described for estimating the fatigue life of 
attitude control thrusters under normal and accelerated loading conditions . 
The method itself, has capabilities which exceed the current thruster 
application. A discussion of possible extensions of the method is 
given in the Concluding Remarks . 

The objective of Volume II of this report is to illustrate In detail the use 
of the RETSCP finite element computer program for engine configurations 
where elastic strains dominate; that is , high cycle fatigue . The 
combustion chamber configuration, operating conditions, and test 
results are described In sufficient detail to conduct the strain analysis 
and compare results. These conditions correspond to an engine which 
was fatigue tested at the NASA Lewis Research Center, The strain 
analysis Is presented In such a manner that the work could be extended 
to other configurations and operating conditions . The fatigue life 
analysis section describes the predictive method employed; and , 
also Indicates methods which could be applied for further studies . 
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ATTITUDE CONTROL THRUSTER CONFIGURATION 
The attitude control thruster configuration to be analyzed is shown in 
Figure 1 from Reference 1. The thruster consists of four major 
components; namely, injector, thrust chamber, propellant valves , 
and igniter assembly. These components are integrated to provide a 
lightweight, compact, high performance thruster, which meets the 
engine operating conditions and the engine design requirements as 
shown in Table I and II respectively. The thruster pulse life design 
requirement is based on 1000 pulses per mission with 100 missions 
per year and a design safety factor of 5 for a total design pulse life 
of 500,000 pulses. 

Although all attitude control thruster components must be capable of 
meeting the above requirements; the combustion thrust chamber, which 
is exposed to a thermal and pressure pulsing environment, deserves the 
most consideration from the standpoint of fatigue. The combustion 
thrust chamber is comprised of injector, chamber, throat section, 
nozzle , torus , liner and flanges which are welded together to form an 


integral unit 



TABLE I 

ENGINE OPERATING CONDITIONS 
(Nominal Design Values) 


Thrust 

Chamber Pressure 
Mixture Ratio 
Nozzle Expansion Ratio 
Propellant Inlet Temperature: 
Hydrogen 
Oxygen 

Propellant Inlet Pressure (to Valve): 
Hydrogen 
Oxygen 

Specific Impulse: 

Steady State 
Pulsing Average 


1500 lb (6672 N) 

300 psla (207 x 10^ N/m^) 

4.0 

40.1 

250°R (139®K) 

375°R (208°K) 

400 psla (276 x lO^ N/m^) 

400 psla (276 x 10^ N/m^) 

435 lb--sec/lb (4266 N-sec/kg) 
I m 

400 lbr“sec/lb (3923 N-sec/kg) 


X 


I 


TABLE II 


ENGINE DESIGN REQUIREMENTS 


Fuel: 


Gaseous hydrogen derived from the 
vaporization of liquid hydrogen. 


Oxidizer: 


Gaseous oxygen derived from the vaporization 
of liquid oxygen. 


Installation: 


Buried within vehicle mold line. 


Maximum External Temperature: 


500°F (260° C) 


Total Life Capability: 


Estimated 2500 minutes. 


Total Number of Firings: 


Estimated 500,000 pulses, plus 25,000 deep 
thermal cycles (full temperature range on 
each component). 


Maximum Single Firing Duration: 


SOO sec. 


Compatibility: 


Reusability: 


Compatible with propellants, test fluids, 
cleaning fluids and environmental contaminants 
for ten-year life requirement. 

To be reusable with minimum servicing and 
refurbishment. 


Service and Maintainability: 


Minimum Impulse Bit (goal): 
Response: 


Reentry Heating (goals): 


Nozzle Scarfing Capability: 


Design for ease of service and maintenance 
when required, 

50 lb-sec (222 N-sec) 

50 mDlisec (time from electrical signal to 
90% thrust). 

30 minutes exposure per mission to the 
following temperatures: 

At nozzle exit plane: 1800°F (982° C) 
At chamber throat: 1200°F (649° C) 

Nozzle skirt shall be easily scarfed beyond 
an expansion ratio of 25:1. 


Weight Go^ for Complete Thruster 15 lb (6804 g) 
Assembly Less Valves: 
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The Injector is a 72 element " premix I” F-O-F triplet design* The face 
of the injector is nade from platelets that contain a fuel '*'1 ■ ‘ cooling 
circuit, the mixing cups, and the 144 "I” shaped fuel massages v.hlch 
discharge approximately 180^ apart at right angles to the oxidizer flow 
(oxidizer flows normal to the face through 72 tubes; fuel flows parallel 
to the face impinging two fuel elements on each oxidizer element). 

The hydrogen flow into the injector comes from the outlet of the chamber 
regenerative passages. 

The thrust chamber incorporates a composite cooling system. The major 
portion of the combustion chamber is regeneratively cooled (80 passages) 
with fuel. The convergent nozzle contains 160 orificed passages which 
direct 21% of the fuel flow to a station 1.50 in. (3.81 cm) upstream of 
the throat where it is introduced as a film coolant. The throat and entire 
divergent no'^zle are of spun Haynes 188 material and operate at the 
adiabatic wall temperature. The gas-side M^all material in the cooled 
area forward of the fuel film coolant injection station is zirconium 
copper alloy which contains rectangular cross-section coolant passages. 
These are closed out by a 22-13-5 stainless steel shell which is brazed 
to the copper liner. The stainless steel shell is Joined to the injector 
and adiabatic wall nozzle by welding . The composite cooling system 
results in the fuel inlet manifold (which is also made from 22-13-5 
material) being located at the entrance to the convergent nozzle . 


The following analysis was conducted specifically at the film cooled 
throat section. The throat section has 1.92 inch (4.88 cm) inside 
diameter with 0.030 inch (0.076 cm) wall thickness, and is fabricated 
from Haynes Alloy No. 188, References 2 and 3. 

The engine was pulse tested at the NASA Lewis Research Center. Tests 
were conducted at 330 psi (228xld^N/mS chamber pressure and the engine 
developed 1100 pounds (4893N) thrust. Total pulse time was 0.15 
seconds of which approximately 0.10 seconds was operation at full 
chamber pressure. The outer wall temperature at the nozzle throat was 
measured, and found to vary from 1040*^F (560°C) during engine firing, 
to 960®F (51^ C) during engine-off; once the uniform engine pulse train 
had been established. Jt was estimated that the temper? .fference 
( A T) across the wall was + 180® F ( +100® C). Obvious ..e inner wall 
is hotter ( + A T) during engine firing. Since the coolant continues to 
flow after shutdown, the inner wall is cooler (-AT) during engine-off. 
These conditions are represented schematically in Figure 2. 

The attitude control thruster was subjected to 50,825 pulses of the above 


type 



Temperature F 

(®c) 


1200 

(648) 


1000 

(538) 


800 

(427) 



Measured Outer- 
Wall Temperature 


Estimated Inner Wall Temperature 
(Assumed 180° Out of Phase) 


Chamber 
Pressure ' psia 

(N/m^) 



Time (seconds) 


Figure 2 . Attitude control thruster test conditions . 
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STRUCTURAL ANALYSIS 

The structural analysis was performed using the RETSCP finite element 
computer program described in Reference 4. It is well known that fatigue 
life depends on the cyclic strain range magnitude through which the 
component is cycled. Reference 5 and 6. Thus, the objective here is 
to predict the stress-strain cycle for the attitude control thruster throat 
station. 

The RETSCP program used the 12-element model of the throat section 
shown in Figure 3. Each element has thickness 0.005 inch (0.013 cm) 
in the z-direction. The model considered is a one degree wedge sector ; 
thus, the section taken is one of 360 symmetric sections. A twelve 
element model was selected to adequately represent the structure. 

Smaller elements were provided near the inner wall, where plastic 
behavior is most likely to occur. 

Boundary conditions in the RETSCP program are in the form of prescribed 
nodal point forces or nodal point displacements . The symmetry of the 
12-elenent model was previously noted. The corresponding boundary 
condition is zero displacement of boundary points normal to the symmetry 
plane, and free movement of nodal points in the radial direction. This 
condition is indicated by the roller notation in Figure 3. The symmetry 
condition along the 1° surface is accommodated by a coordinate trans- 
formation within the RETSCP program. 
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Figure 3, Finite element model 




The boundary conditions in the z-direction require special consideration. 

A very thin section was taken at the throat, so that the rocket engine 
contour curvature can be neglected. That is, constant thickness elements 
are used. Under load, the section remains plane; and also, the total 
load on the section is equal to the engine thrust load. This corresponds 
to the generalized plane strain condition. 


Under these conditions an approximation for the axial boundary condition 
is obtained from the generalized plane strain equation of elasticity, 
namely (Reference 7) 

I r T 


®2 * E ®z -V ( Ojj + Oy ) (1) 

The stress in the x-direction is computed from the hoop pressure load. 
Stresses in the y-direction are negligible. The z-component of stress is 
determined by the thrust load. Thus, the nodal points along the face z 
* 0 (odd numbers) are fixed and the nodal points on the opposite face 
(even numbers) are given the prescribed displacement determined from 
Equation (1). Since the nodal point forces in the z-direction are printed 
output for RETSCP, a post run evaluation of the net z force is made. If 
the result differs substantially from the thrust force value; then the 
z-direction displacement is adjusted and an iteration is performed. The 
reference temperature is taken to be the mean temperature; thus , the axial 
strain computed using equation (1) need not be adjusted for thermal growth. 
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4 2 

During the engine firing « the throat pressure load is 175 psi (121x10 N/m ). 
The mean wall temperature is 1000°F (538°C) and a linear temperature 
gradient is assumed. That is, the outer wall is at 910^F (488^0 and the 
inner wall at 1090°F (588°C) with constant gradient through the thickness. 

The material properties which were used in the calculation are listed in 
Table III, Reference 2 and 3. Input and output data for the hot engine 
firing pulse are given in Appexdixes B and C. See Reference 4 for card 
content and data format. 

The important point to be noted , from the output data in Appendix C , is 
that the entire structure behaves elastically. Thus, there will be no 
residual strain influence on the engine-off portion of the stress-strain 
cycle . 

The hoop stress component distribution is plotted in Figure 4. The result 
is nearly linear and corresponds to a linear temperature gradient across 
a thin wall cylinder. The hot inner wall corresponding to the engine run 
pulse is in compression; whereas, the cooler outer wall is in tension. 

The stress result for this engine run pulse is non-symmetrical with respect 
to the mean radius due to the chamber pressure load. The structural 
behavior is dominated by the thermal loading . 
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TABLE III 


HAYNES 188 MATERIAL PROPERTIES DATA 


Property 

English Units 
(RETSCP) 

Metric 

Units 

Thermal Expansion Coefficient: 

8.4x10"® in/ln-°F 

15.1x10 ®m/m-®C 

Modulus of Elasticity: 

27.6x10® psi 

10 2 
19.0x10 N/m 

Poisson's Ratio: 

0.30 

0.30 

Plastic Modulus Ratio: 

0.10 

0.10 

Yield Stress * : 



Reference Yield ( a^) 

43,800 psi 

8 2 
3.02x10 N/m 

Yield Slope ( ) 

5.0 psi/ °F 

+4 , 2 o 

6.2x10 N/m -°C 


* . Oq - X, 6r 


Where: 0 = T - T 
r r 

Tj. = 1000°F (538°C) 
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Hood Stress - psi 
(10®x N/m^) 


40,000 

(2.758) 


Engine-Run 

Pulse 


P.nalne-Off 




The engine-off portion of the thruster pulse was computed in a manner 
analogous to that described above. That is, a linear temperature gradient 
between the opposite limits of 1090°F (588°C) on the outer wall and 
910°F (488°C) on the inner wall. Also, the throat station chamber pressure 
is now zero and there is zero net thrust. The appropriate boundary condition, 
in the z-direction, is zero strain. 


Results for the engine-off case are similar to the run condition; except 
that now the outer wall is in compression and the inner wall in tension. 
Again, the entire structure remains elastic. The hoop stress component 
distribution is plotted in Figure 4. Note, that in the absence of pressure 
loading, the stress distribution is nearly symmetrical with respect to the 
mean radius. The stress magnitudes are comparable for the engine-run 
and engine-off pulses. It should be noted that the outer wall stress and 
the inner wall stress alternates between tension and compression as the 
engine is cycled. 


The effective, or total, strain is related to the effective stress via the 
material stress-strain data. Effective stress is defined as follows 
(Reference 4): 


" k 


1/2 


( 2 ) 
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The distribution of effective strain is plotted in Figure 5 for both engine - 
run and engine-off pulses. Since the calculated effective stress can be 
positive or negative, it is necessary to consider the signs of the average 
X - and z - stress components to determine if the effective strain in 
Figure 5 is positive or negative (the y-stress is negligible when compared 
to the X - or z - stress components). The effective strain is considered 
positive when both the average x - and z - stress components are positive. 
Likewise, the effective stress is considered negative when both the 
average x - stress and z - stress components are negative. Near the 
zero strain axis, the x - and z - stress component signs differ so that 
the above criteria does not apply. In that case, smooth curves were 
drawn through the plotted points which met the sign criteria on either side 
of the axis. Associated with the tension-compression stress reversal, 
there is a corresponding strain reversal. Thus, a strain range is estab- 
lished for the engine pulse mode . 






Effective 

Strain 

(%) 

15 

Engine-Off 

Pulse 



i 

! 

i 

? 

I 

1 


Figure 5. Effective strain d 
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FATIGUE ANALYSIS 


* 


Prediction of chamber fatigue life Is made by using the computed strain 
results In conjunction with uniaxial fatigue data. 


In the case of Haynes 188, suitable fatigue test data are not available. 
Thus, strain-life curves are generated by the Universal Slopes method. 
Reference 6; namely. 


et 


3.5 - -0,12 , 

— F — N, + In 



0.60 _ -0,60 


f} 


f 


(3) 


where, the ultimate tensile strength , modulus of elasticity E, and 
reduction In area RA are based on the short term tensile test. 


Now, for elevated test temperatures, the Universal Slopes prediction Is 
modified. Reference 8. First, equation (3) Is based on the short term 
properties at the elevated temperature. Then, to account for reduced life 
at elevated temperature, the “minimum" life N| Is defined as one tenth of 
the value In equation (3), N^ /lO. Average life Is twice this minimum 
value. This procedure Is based on a comparison with a great deal of 
experimental data. 

Minimum life values were computed per the above m<»thod for temperatures 
In the range 700 - 1300°F (371-704° C), Results are presented In Figure 6. 
The curves coincide over the range considered. Typically at 1000°F (538®C), 
the fatigue life curve Is based on E ■ 27.6xl0^psl (19.0xl0^^N/m^), 

8 9 

Oy ■ 107,200 pfi (7.39x10 N/m^) , and RA ■ 0.60. 
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Strain Range 



Figure 6. Minimum fatigue life for Haynes 188 at elevated temperature. 





4 


The linear cummulative damage law« Reference 9, is useful when the 
loading produces varying strain range values. According to that law« 
if nj and nj^ cycles are applied at conditions for which the fatigue life 
would be Nr and N. respectively; then, the damage is accumulated 
linearly as follows: 


N, 


N, 


(4) 


Failure occurs when the cumulated damage totals unity. 


For the case at hand, mean or typical loading conditions are assumed. 
Thus , the fatigue life is obtained directly from the material fatigue 
test data. Had the strain calculation been performed for various loading 
conditions; the life would be computed using equation (4). Much work 
has been done in the area of fatigue damage , along both the experimental 
and analytical lines . The topic is touched on briefly here only to Indicate 
how an accurate loading spectrum would be treated. 


Now , for the conditions at hand; fatigue life is determined directly by 
entering Figure 6 with the maximum value of strain range . The resulting 
predicted minimum life is 220,000 cycles. Recall that, the number of 
pulses, which the nozzle throat was subjected to without failure, 
was 50,82s. 


Ifc 1 1 V'^tuSton 
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CONCLUDING REMARKS 


The detailed life analysis presented herein was conducted for configuration 
and operating conditions which coirespond directly to actual Integrated 
thruster assembly test conditions. Thus, comparison of predicted and 
observed fatigue behavior is possible. The computed strain range was 
used to predict the minimum number of cycles to failure. The predicted 
minimum life was in fact greater than the number of pulses to which the 
thruster was subjected without failure. 

The RETF P program has three-dimensional and elasto-plastic capability. 
Thus , the program and analytical procedure have potential for extended 
applications in attitude control thruster evaluation. For example, complex 
configurations of other locations within the engine could be analyzed in 
detail. These include the regeneratively cooled combustion chamber, 
and the film coolant exit lip. 

Parametric studies, at the nozzle throat, could also be conducted. That 
is, the variation in fatigue life with chamber pressure, O/F ratio, and 
wall thickness could be determined. Also, by adding several rows of 
elements , the effects of axial temperature gradient and nozzle throat 
curvature could be studied. 

Application of the above techniques would not only generate more efficient 
thruster designs; but would also increase understanding of the parameters 
which affect engine fatigue life. 
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APPENDIX A SYMBOLS 

D Fdtlgue damage fraction 

E Modulus of elasticity 

n Number of loading cycles 

N£ Minimum fatigue life 

Universal slopes cycles to failure 
RA Reduction Ir area 

T Temperature 

Reference temperature 
x,y,z Cartesian coordinates 


c Strain 

0 Temperature difference , T -T^. 

o Stress 



0 Effective stress 

e 

OxtOy,oi Normal stress components 
Oy Yield stress 

0^ Yield stress at reference temperature 

Cy Ultimate tensile strength 

X| Yield slope with temperature 

j Shear stress 

V Poisson's ratio 

^ 25 - 




APPENDIX B SAMPLE RETSCP INPUT DATA 


The following .nput listings are from the engine-run portion of the 
test cycle* The data are Included as programming guidlines. For 
complete description of input data and format see Reference 4. 
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APPENDIX C SAMPLE RETSCP OUTPUT DATA 


The following output data are from the engine-run portion of the test 
cycle. The nodal point displacements (inches), stress components (psi), 
and the yield check summary (Inches/inch or psi where applicable) are 
Included. For complete descirption of output data and formats see 
Reference 4. 
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